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ABSTRACT 

We report on the results of an analysis of Chandra, XMM-Newton and new GMRT data of the X-ray bright 
compact group of galaxies HCG 62, which is one of the few groups known to possess clear, small X-ray cavities 
in the inner regions. This is part of an ongoing X-ray /low-frequency radio study of 18 groups, initially chosen 
for the availability of good-quality X-ray data and evidence for AGN/hot gas interaction. At higher frequency 
(1.4 GHz), the HCG 62 cavity system shows minimal if any radio emission, but the new GM/?r observations 
at 235 MHz and 610 MHz clearly detect extended low-frequency emission from radio lobes corresponding to 
the cavities. By means of the synergy of X-ray and low-frequency radio observations, we compare and discuss 
the morphology, luminosity and pressure of the gas and of the radio source. We find that the radio source is 
radiatively inefficient, with a ratio of radio luminosity to mechanical cavity power of ^ lO"'*, and that the radio 
pressure of the lobes is about one order of magnitude lower than the X-ray pressure of the surrounding thermal 
gas. Thanks to the high spatial resolution of the Chandra surface brightness and temperature profiles, we also 
identify a shock front located at 36 kpc to the south-west of the group center, close to the southern radio lobe, 
with a Mach number L5 and a total power which is about one order of magnitude higher than the cavity 
power. Such a shock may have heated the gas in the southern region, as indicated by the temperature map. 
The shock may also explain the arc-like region of enriched gas seen in the iron abundance map, as this may be 
produced by a non-Maxwellian electron distribution near its front. 

Subject headings: galaxies:clusters:general - galaxies: clusters: individual (HCG 62) - cooling flows - inter- 
galactic medium - galaxies:active - X-rays: galaxies: clusters 



1. INTRODUCTION 

Many possibilities have been proposed in the last decade to 
solve the so-called "cooling flow problem" in the hot atmo- 
spheres of galaxy clusters and groups., i.e., the lack of evi- 
dence for central gas cooling to very low temperatures at the 
rates predicted (see, e.g., the review by Peterson & Fabian 
2006). Among these, feedback by the central active galac- 
tic nucleus (AGN) appears to be the most promising solution. 
Identifying radio galaxies (galaxies whose AGN produce the 
radio-emitting jets) as a primary source of feedback in clusters 
has been one of the major achievements of the current gener- 
ation of X-ray observatories. Although gas does cool at least 
through part of the X-ray emitting temperature band (10^-10** 
K) in galaxy clusters, Chandra and XMM-Newton have shown 
that there is not a significant amount of gas cooling below 
about one third of its original temperature (e.g., Peterson et 
al. 2003; Kaastra et al. 2004). These observations also re- 
veal highly disturbed structures in the cores of many clusters, 
including shocks, cavities and sharp density discontinuities. 
At radio wavelengths, it is clear that AGN jets are the cause 
of many of these disturbances. The most typical configura- 
tion is for jets from the central dominant elliptical of a cluster 
to extend outwards and inflate lobes of radio-emitting plasma. 
These lobes push aside the X-ray emitting gas of the cluster at- 
mosphere, thus leaving apparent cavities in the X-ray images 
(e.g., for a review: McNamara & Nulsen 2007 and references 
therein). 

There are already several in-depth analyses of individual 
rich clusters (e.g., Perseus: Fabian et al. 2006; Centaurus: 
Sanders et al. 2008; MS0735.6H-7421: Gitti et al. 2007; 



A2052: Blanton et al. 2009) and studies of cluster samples 
(Birzan et al. 2004, 2008; Dunn et al. 2005; Dunn & Fabian 
2006; Rafferty et al. 2006; Diehl et al. 2008). However, 
the differences between groups and clusters imply that the ex- 
isting studies on feedback tell us little about how it operates 
in groups. Furthermore, due to the shallower group poten- 
tial, the AGN outburst is expected to have a large impact on 
the intragroup medium. It is therefore essential to study indi- 
vidual groups and group samples in order to understand how 
feedback has influenced the thermal history of galaxies and 
the intragroup medium, and thus of most of the baryons in the 
Universe. Such investigations have been undertaken only re- 
cently (e.g., Johnson et al. 2009; Sun 2009; Gastaldello 2008, 
2009; Giacintucci et al. in prep.; O' Sullivan et al. in prep.). 

The work presented here is part of an ongoing project aimed 
at combining X-ray and low-frequency radio observations of 
galaxy groups. In particular, we have selected a sample of 
18 galaxy groups based on the presence of signs of inter- 
action between the hot gas and the central AGN. For these 
groups, which all have good quality X-ray data in the archives 
of Chandra and/or XMM-Newton, we have obtained new ra- 
dio data at the Giant Metrewave Radio Telescope (GMRT) at 
235 MHz, 327 MHz, and 610 MHz (Giacintucci et al. 2008, 
2009; Raychaudhury et al. 2009, Giacintucci et al. in prep.). 

HCG 62 (z=0.0137) is the X-ray brightest of the Hick- 
son compact groups and has been extensively studied with 
ROSAT, ASCA, Chandra, XMM-Newton and Suzaku (Ponman 
& Bertram 1993; Finoguenov & Ponman 1999; Buote 2000; 
Morita et al. 2006; Gu et al. 2007; Tokoi et al. 2008; Sanders 
et al. 2009). The central region of this group is dominated 
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by four early-type galaxies, two of which are likely interact- 
ing as indicated by optical studies (Valluri & Anupama 1996; 
Spavone et al. 2006). This was the first galaxy group with a 
clear detection of inner cavities (Vrtilek et al. 2002). The ex- 
isting 1 .4 GHz VLA observations mainly show the emission 
from the compact radio source and indicate only some hints of 
extended radio emission toward the southern cavity (Vrtilek et 
al. 2002). Owing to the poor radio images then available, the 
HCG 62 cavity system was classified as a "radio ghost" in the 
sample of Birzan et al. (2004). 

We present here new low frequency GMRT radio observa- 
tions of HCG 62 that, together with the existing X-ray Chan- 
dra and XMM-Newton observations, allow us to study the X- 
ray /radio interaction and investigate the AGN feedback in this 
group. This paper is organized as follows. In ^we describe 
the data sets used for our investigation and the data reduction 
process, in Owe show the X-ray and radio morphologies of 
the group inner region, and in ^we present the X-ray proper- 
ties derived from the spectral analysis. In ^5.1| we investigate 
the interaction of the radio plasma with the hot gas in terms 
of energetics 05.1.1I ) and pressure balance 05.1.2I ). In ^5.21 
we analyze the surface brightness discontinuities, finding ev- 
idence for the detection of weak shocks. Finally, in ^5.3| we 
present the 2-D distribution of temperature ( ^5.3. Il l and iron 
abundance 05.3.2b . The summary of our main results is given 
in|6] 

With //() = 70 km s"' Mpc"', and 9.m=1-^k = 0.3, the 
luminosity distance to HCG 62 is 59 Mpc and 1 arcsec corre- 
sponds to 0.28 kpc in the rest frame of the group. The radio 
spectral index a is defined such as Si, oc where Si, is the 
flux density at the frequency f. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

HCG 62 was observed with Chandra for ^ 48.5 ks on 2000 
January 25 (ObsID 921) with ACIS-S in imaging mode op- 
erating at the focal plane temperature of -110°C. The data 
were reprocessed with CIAO 4.1 using CALDB 4.1.0 and 
corrected for known time-dependent gain problems following 
techniques similar to those described in the Chandra analysis 
threads'. No charge transfer inefficiency (CTI) correction is 
available for data taken at -1 10°C. Screening of the event files 
was applied to filter out strong background flares. Blank-sky 
background files, filtered in the same manner as in the HCG 
62 image and normalized to the count rate of the source image 
in the 10.0-12.0 keV band, were used for background subtrac- 
tion. The final exposure time is 46.5 ks. We only use data 
from the S3 CCD since it covers the central part of the group 
emission where the cavity system and the radio source are lo- 
cated. We identified and removed the point sources on S3 us- 
ing the CIAO task WAVDETECT, with the detection threshold 
set to the default value of lO"*". 

Updated gain correction at -110°C only become available 
with CIAO 4.0 CALDB version, so previous work has 
been limited by calibration uncertainties. In particular, 
spectra were extracted using die SPECEXTRACT task and, 
as recommended for -110°C data^, new RMF files were 
computed using the task MKACISRMF with the latest gain 
file, acisD2000-01-29p2_respN0005 . fits. 



' http://cxc.harvard.edu/ciao/threads/index.htmI 
^ |http://cxc.harvard.edu/ciao/why/l lO.html | 



2.2. XMM-Newton 

HCG 62 was observed by XMM-Newton in June 2007 dur- 
ing revolution 1382 (ObsID 0504780501, hereafter #501) and 
1383 (ObsID 0504780601, hereafter #601) for nominal expo- 
sure times of 122.5 ks and 32 ks, respectively. In this pa- 
per, only data from the European Photon Imaging Camera 
(EPIC) are presented. The MOS detectors were operated in 
Full Frame Mode and the PN detector in Extended Full Frame 
Mode, both with the MEDIUM filter. We used the SASv8.0.0 
processing tasks emchain and epchain to generate cali- 
brated event files from raw data. Throughout this analysis 
single and double pixel events for the PN data (PATTERN 
<=4) were selected, while for the MOS data sets the PAT- 
TERNS 0-12 were used. The removal of bright pixels and hot 
columns was done by applying the expression (FLAG==0). 
To reject soft proton flares, we generated a light curve in the 
10.0-12.0 keV band where the emission is dominated by the 
particle-induced background, and excluded all the intervals 
of exposure time having a count rate higher than a certain 
threshold value (the chosen threshold values are 0.20 (0.18) 
cps for MOS and 0.35 (0.32) cps for PN in obs. #501 (#601)). 
The total remaining exposure times after cleaning are 112.0 
ks for MOSl, 97.6 ks for MOS2 and 60.2 ks for PN. Starting 
from the output of the SAS detection source task, we made 
a visual selection on a wide energy band MOS & PN image 
of point sources in the field-of-view. Events from these re- 
gions were excluded directly from each event list. The source 
and background events were corrected for vignetting using 
the weighted method described in Arnaud et al. (2001), the 
weight coefficients being tabulated in the event list with the 
SAS task evigweight. This allows us to use the on-axis 
response matrices and effective areas, computed in the central 
30" with the tasks rmf gen and arf gen. 

2.2.1. Background treatment 

The background estimates are derived using blank-sky ob- 
servations consisting of a superposition of pointed observa- 
tions that have been processed with SAS version 7.1.0 (Carter 
and Read 2007). In particular, for each camera, we obtained 
a blank-sky template tailored to the characteristics of HCG 
62 specific observation by submitting a XMM-Newton EPIC 
Background Blank Sky Products Request Form^^ with the re- 
quest"* of a Galactic column in the range 2 to 5 x 10^"cm"^. 

The blank-sky background events were then selected using 
the same selection criteria (such as PATTERN, FLAG, etc.), 
intensity filter (for flare rejection) and point source removal 
used for the observation events. This yields final exposure 
times for the blank fields of 1.0 Ms for MOSl, 1.6 Ms for 
MOS2 and 99.3 ks for PN. Since the cosmic ray induced 
background might change slightly with time, we computed 
the ratio of the total count rates in the high energy band (10.0- 
12.0 keV). The obtained normalization factors (1.258 (1.203), 
1.316 (1.246), 1.353 (1.175) for MOSl, MOS2 and PN, re- 
spectively, in obs. #501 (#601)) were then used to renormalize 
the blank field data. The blank-sky background files were fi- 
nally recast in order to have the same sky coordinates as HCG 
62. 

The usual approach to perform the background subtraction 

' http://www.star.le.ac.uk/~jac48/BG/UserRequest/blankskyform.html 
The selection criteria adopted are: Filter: Medium, Model: Full-Frame 
(Extended Full-Frame for PN), Type: Ghosted, Galactic column: 2.0 — 5.0 X 
lO^^cm"^, Equatorial coordinates: (193.274,-9.2036) degrees, radius = 180 
degrees. 



Fig. 1 . — Top: Raw (left) and unshaq3 masked (right) 0.5-2.0 keV ACIS-S image of the central region of HCG 62. Bottom: Raw (left) and unsharp masked 
(right) 0.5-2.0 keV MOS1-I-MOS2-I-PN mosaic image of the central region of HCG 62. The unsharp masked images are produced b y su btracting a large-scale (20 
pixels) smoothed image from a small-scale (5 pixels) smoothed image. The arrows indicate the features discussed in the text (see j|3.U . In particular, the dashed 
arc in top panels indicates the position of the SW shock, which is non visible clearly in the images. In all panels, the box is 6' X 6' (101 kpc by 101 kpc). North 
is up and east is left. 



is described in full detail in Arnaud et al. (2002). This pro- 
cedure consists of two steps. In the first step, for each spec- 
trum extracted from the observation event list an equivalent 
spectrum is extracted from the corresponding blank-field file 
and then subtracted from it. This allows us to remove the 
particle background. However, if the background in the ob- 
servation region is different from the average background in 
blank field data, this step could leave a residual background 
component. The residual background component is estimated 
by using blank-field-subtracted data in a region supposedly 
free of source emission (in particular we considered an an- 
nulus lying between 13 and 14 arcmin) and then subtracted 
in a second step from each MOS and PN spectrum. We note 
that while some group emission might still fall in the annulus 



chosen to estimate the residual background component, our 
spectral analysis is not affected as the central region we are 
interested in is very bright and source dominated. 

2.3. GMRT 

HCG 62 was observed by the GMRT in February 2008 at 
frequencies of 610 MHz and 235 MHz (project 13SGa01) for 
a total of 125 and 250 minutes, respectively. The data were 
acquired in spectral line mode with 128 channels/band and a 
bandwidth of 16 MHz at 610 MHz and 8 MHz at 235 MHz. 
Data reduction was done using the NRAO AIPS (Astronom- 
ical Image Processing System) package. Accurate editing of 
the uv data was applied to identify and remove bad data. Af- 
ter bandpass calibration, the central channels were averaged 



4 





Fig. 2. — Left: GMRTbiO MHz contours overlaid on the unsharp masked 0.5-2.0 keV Chandra image (same as in top riglit of Fig. [T). Tlie beam size is 14" 
by 14" and the lowest contour is shown at 3cr = 0.15 mJy/beam. Right: GMRT 235 MHz contours overlaid on the smoothed 0.5-2.0 keV Chandra image. The 
beam size is 14.4" by 13" and the r.m.s. noise is 170/iJy/beam. The contours are spaced by a factor 2 stalling from the lowest level of 0.45 mJy/beam. In both 
panels, the cross indicates the position of the radio core. North is up and east is left. 

to 6 channels of - 2 MHz at 610 MHz and - 1 MHz at 235 
MHz. After further careful editing in the averaged datasets, 
images were produced using the standard procedure (calibra- 
tion, Fourier- Transform, Clean and Restore) and the wide- 
field imaging technique. Phase-only self-calibration was ap- 
plied to remove residual phase variations. The rms noise level 
(Icr) in the final images is 50 /iJy/beam at 610 MHz and 170 
/ijy/beam at 235 MHz. Following Chandra et al. (2004), we 
estimated that the amplitude calibration error at 610 MHz is 
^5%, while a calibration error of 8% was assumed at 235 
MHz. The calibration uncertainties were included in the flux 
density and spectral index errors determined in 33.21 The 
GMRT observations of HCG 62 will be presented and dis- 
cussed in more detail in Giacintucci et al. (in prep.). 

3. THE DISTURBED INNER REGION OF HCG 62 

3.1. X-ray morphology 

The raw 0.5-2.0 keV ACIS-S image in Fig. [T] (top left 
panel) shows that the hot gas in the central region of HCG 
62 is not smoothly distributed but is instead perturbed with 
cavities and edges. In particular, as best revealed by the un- 
sharp masked image shown in Fig. [T] (top right panel), we 
confirm the presence of two clear cavities, first discovered by 
Vrtilek et al. (2002), and we also notice the presence of sur- 
face brightness discontinuities or "fronts" which suggest the 
presence of large-scale gas motions or shocks in the group 
halo. In particular, by means of the surface brightness pro- 
files presented in 35.21 we identify a front at 36 kpc 2') 
from the center toward the south-west (SW) at a distance of 
about 10 kpc 35") from the southern radio lobe, and one 
at 20 kpc 70") from the center toward the north-east (NE) 
close to the N cavity. These features are significant at the 5a 
and 15(7 level, respectively. In the images there is also a hint 
of an edge about the same distance SW of the group center 
as the front to the NE. The spectral analysis presented be- 



low indicates that the outer SW edge is a shock front ( 35.21 ). 
As already estimated by previous studies (Birzan et al. 2004, 
Rafferty et al. 2006, Morita et al. 2006) the two cavities have 
similar sizes, with diameters of 10 kpc, and are located at 
a projected distance of ~ 8.5 kpc from the group center. We 
extracted the surface brightness profiles along and orthogonal 
to the cavities, and at the radius of 8.5 kpc we estimate a 
brightness decrement of ~ 30% along the cavities relative to 
the orthogonal direction. 

Because of the poorer XMM-Newton spatial resolution, the 
cavities are not evident in the raw 0.5-2.0 keV image shown 
in Fig. [T](bottom left panel). However, they do appear clearly 
in the unsharp masked image (Fig. [T] bottom right panel), at 
positions coincident with those determined by the Chandra 
data. 

3.2. Radio properties 

The existing higher frequency (1.4 GHz) observations of 
HCG 62 show the central compact radio source with some 
hint of emission extending toward the southern cavity (Vr- 
tilek et al. 2002). With the new GMRT observations at 235 
MHz and 610 MHz we detect clearly extended emission em- 
anating from the core in the typical form of a bipolar flow. 
The two radio lobes point toward the N and S directions, fill- 
ing the cavities. The overlay of the 610 MHz radio contours 
on the unsharp masked Chandra image is shown in Fig. |2] 
(left), whereas the overlay of the 235 MHz radio contours on 
the smoothed Chandra image is shown in Fig. |2] (right). It 
is evident that far more extensive structures become visible at 
lower frequencies: the radio emission at 235 MHz is more ex- 
tended with two faint regions located outside the cavities and 
apparently bent toward the east (E). In particular, we identify 
inner lobes clearly filling the well defined X-ray cavities (best 
visible in left panel of Fig. |2|, and outer lobes having no as- 
sociated X-ray cavities (best visible in right panel of Fig. |2]). 
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TABLE 1 

XMM-Newton PROJECTED SPECTRAL Analysis 



TABLE 2 

XMM-Newton Deprojected SPECTRAL Analysis 
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Note. — Results of the XMM-Newton spectral fitting in concentric 360° an- 
nular regions in the 0.7-3.0 keV energy range using the absorbed vapec ther- 
mal plasma model and fixing the absorbing column density to the Galactic value 
(A'h = 3.31 X 10^0 cni-^) xhe temperature (in keV) and abundances (in fraction of 
the solar value, Anders & Grevesse 1989) are left as free parameters. Error bars 
are at the 90% confidence levels on a single parameter of interest. The first column 
indicate the number region. The delimiting radii are: reg.#l: 0-15 kpc (0-55"), 
reg.#2: 15-28 kpc (55"-100"), reg.#3: 28-42 kpc (100"-150"), reg.#4: 42-56 kpc 
(150" -200"), reg.#5: 56-70 kpc (200"-250"), reg.#6: 70-126 kpc (250"-450"). 

The connection of the radio source with the X-ray features 
will be discussed in ^ 

At 610 MHz, the source has a total flux density of 13.5 ±0.7 
mJy. At 235 MHz, the total flux density is 42.4 ±3.4 mJy. 
Neglecting the K-correction term (for 0.3 < a < 1.7 the K- 
correction is < 1%, Petrosian & Dickey 1973), the monochro- 
matic radio power at each frequency is calculated as 



Note. — Results of the XMM-Newton deprojection analysis in con- 
cenfi'ic 360° annular regions in the 0.7-3.0 keV energy range using the 
XSPEC pro jet Xwabs Xmekal model and fixing the absorbing col- 
umn density to the Galactic value {Nn = 3.31 X 10^" cm"-). The fit gives 
X^/dof = 19848/4463. The first column indicate the number region. The 
delimiting radii are the same as in Table[T] 



(1) 



where Dl is the luminosity distance, which gives PeiOMHz = 
(5.7±0.3) X lO^i W Hz-' and P235MHZ = (1.8±0.1) x 10^2 
W Hz"'. This radio source is thus classified as a weak FR-I, 
and its spectral index 023*5 mhz ~ 1.2 ±0.1 is steep compared 
to that of typical radio galaxies. 

Since the outer radio lobes are only detected at one fre- 
quency in our GM/?r observations, we are unable to perform 
a detailed radio spectral analysis. New, deeper GMRT obser- 
vations at 235, 327 and 610 MHz, scheduled during Cycle 17, 
will allow us to study with unprecedented detail the spectral 
properties of the source. Based on the analysis of the inte- 
grated radio spectrum and spectral index distribution, an es- 
timate of the radiative ages of the electron population can be 
obtained (e.g., Giacintucci et al. 2007, 2008 and references 
therein) for both pairs of radio lobes. This will help to shed 
light on the relationship of the outer lobes to the inner ones. 

4. GLOBAL X-RAY PROPERTIES 

Throughout the analysis, a single spectrum for each instru- 
ment and each dataset was extracted for each region of inter- 
est. The relative normalizations of the MOS and PN spectra 
were left free when fitted simultaneously. Spectral fitting was 
performed in XSPEC version 12.3.1. Abundances were mea- 
sured relative to the abundance ratios of Anders & Grevesse 
(1989). A galactic hydrogen column of 3.31 x lO^^cm'^ 
(Dickey & Lockman 1990) and a redshift of 0.0137 was as- 
sumed and, unless otherwise stated, the reported errors are at 
the 90% confidence level. Spectra were grouped to 25 counts 
per bin and the energy range 0.7-3.0 keV was adopted. 

We first tested in detail the consistency between Chandra 
and the three XMM-Newton cameras by extracting a global 
spectrum from all events lying in the annulus 1-3 arcmin from 
the group emission peak and fitting these spectra separately 
with an absorbed apec model. We find good agreement be- 
tween the three XMM-Newton cameras {kT = 1 .25;|:o [Jj 



for MOSl, 1.21!!]:!]^ keV for MOS2, l.2?,^l\ keV for PN, 
with a value of 1 .23^q q[ for the combined fit), whereas there 
is a discrepancy of about 0. 1 keV with the Chandra measure- 
ment {kT = 1 .32^j}Q[ keV). However, this is consistent with the 
systematic cross-calibration uncertainties (David 2009). Due 
to the systematic offset, we analyzed Chandra and XMM- 
Newton data separately. In particular, we present here the re- 
sults obtained by the analysis of XMM-Newton spectra only 
because they comprise most of the counts. 

4.1. Azimuthally Averaged Projected Temperature and 
Metallicity Profiles 

We produced projected radial temperature and metallicity 
profiles by extracting spectra in six circular annuli centered 
on the peak of the X-ray emission. The annular regions are 
described in Table [T] The data from the three XMM-Newton 
cameras were fitted simultaneously to an absorbed vapec 
thermal plasma model, where O, Si, S, and Fe are treated as 
free parameters with the Ni abundance linked to the Fe abun- 
dance and the Mg and Ne abundances linked to the O abun- 
dance. The best-fitting parameter values and 90% confidence 
levels derived from the fits to the 360° annular spectra are 
summarized in Table [T] 

The azimuthally averaged projected temperature profile 
(Fig. [3] left) shows a positive temperature gradient at small 
radii, with temperatures increasing up to ^ 1 .5 keV at a radius 
of ^ 50 kpc. This behavior is similar to that observed in other 
groups (e.g., Gastaldello et al. 2007; Finoguenov et al. 2007; 
Rasmussen & Ponman 2007; Sun et al. 2009; O'Suflivan et 
al. in prep.), although we do not observe the typical decline 
in temperature at larger radii. Previous ROSAT observations 
of HCG 62 indicate a temperature decline beyond ^120 kpc 
(Ponman & Bertram 1993), so our result is consistent with 
earlier literature. 

The azimuthally averaged projected abundance profiles 
(Fig. [3] right) do not show strong gradients, but are instead 
relatively flat around values of ~0.5 solar for both Si and S, 
and '--^0.3 solar for Fe. We were able to constrain O only in 
the central region, with a value of ^0.4 solar. 

4.2. Deprojection Analysis 

To correct for the effect of projection along the Une of sight, 
we also performed a deprojection analysis on the same an- 
nular spectra used in ^4.1l bv adopting the XSPEC pro jet 
model. Under the assumption of ellipsoidal (in our specific 
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Fig. 3. — Left: Azimuthally averaged projected temperature profile measured witli XMM-Newton and reduced = X^/dof of the fit. Right: Azimuthally 
averaged projected abundance profiles of Si, S, and Fe measured with XMM— Newton. 
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Fig. 4. — Left: Azimuthally averaged deprojected profiles of temperature (top) and density (bottom) measured with XMM— Newton. Right: Azimuthally 
averaged deprojected profiles of pressure (top) and entropy (bottom) measured with XMM~Newton. 



case, spherical) shells of emission, this model calculates the 
geometric weighting factor, according to which the emission 
is redistributed amongst the projected annuli. As expected, 
the deprojected central temperature is lower than the pro- 
jected one, since in the projected fits the spectrum of the cen- 
tral annulus is contaminated by hotter emission along the line 
of sight. We also note that the low value that we measure 
probably represents the very coolest gas in the core, owing 
to emission weighting in the large bin adopted to account for 
the XMM PSF. Due to the limited photon statistics relative 
to the high number of parameters in the pro jet model we 
could only poorly constrain the deprojected abundance pro- 
files, which therefore are not shown. 



We also estimate various quantities derived from the depro- 
jected spectral fits. The electron density rig is obtained from 
the Emission Integral EI = J tif-n^dV given by the vapec nor- 
malization: 10"''*£'//(47r[DA(l+z)]^). We assume Hp = 0.82ne 
in the ionized intra-cluster plasma. By starting from the de- 
projected density and temperature values, we can then calcu- 
late the gas pressure as p = nkT, where we assume n = In^, 
and the gas entropy from the commonly adopted definition 
S = kT nj'^^. The results are reported in Table|2]and the corre- 
sponding deprojected temperature, density, pressure and en- 
tropy profiles are shown in Fig. |4] 
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Fig. 5. — Top, left: Overlaid on the smoothed 0.5-2.0 keV XMM~Newton mosaic image are the sectors used for the spectral analysis (see j|4.3l for details). 610 
MHz GMSrradio contours are also shown to locate the position of the inner radio lobes. Top, right: Projected temperature (top panel) and Fe abundance (bottom 
panel) profiles measured with XMM-Newton along sector NE (red circles), sector SW (red squares), sector WN (black triangles) and sector ES (black pentagons). 
In red and black are indicated the profiles measured along and orthogonal to the cavities, respectively. Bottom: "Wedge" temperature (left) and Fe abundance 
(right) maps derived from the XMM-Newton projected spectral fits in the sectoral annuli detailed in top, left figure. For comparison, Chandra contours and 610 
MHz GMRT contours are overlaid in black and green, respectively. 

4.3. Azimuthal Variations in Projected Temperature and 
Metallicity Profiles 

The interaction of the radio source with the thermal gas 
can produce significant effects on the X-ray properties of the 
group halo, that are thus expected to show azimuthal varia- 
tions. To investigate this, we performed a spectral analysis in 
annular sectors each having angular aperture of 45°. Start- 
ing from the W direction with position angle (P.A.) = 0° and 
counting counterclockwise, the sectors are labeled as: WN 
(west-north), NW (north-west), NE (north-east), EN (east- 
north), ES (east-south), SE (south-east), SW (south-west). 



WS (west-south). Each sector was divided in the same six 
annular regions used in ^4.1[ obtaining the regions shown in 
Fig. |5](left panel). We produced projected radial temperature 
and metallicity profiles by extracting spectra in these regions. 
In particular, we compared the radial properties along and or- 
thogonal to the radio lobes and cavity system. The northern 
and southern cavities fall entirely within the first annulus in 
the NE and SW sectors, respectively. The profiles derived 
along the different sectors are shown in Fig. |5] (top right), 
while the temperature and metallicity wedge maps are shown 
in Fig. |5] (bottom left and right, respectively). The gas prop- 
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erties do not appear to show variations obviously associated 
with the cavities and radio lobes. This indicates that the spec- 
tral analysis in large sectors is probably not the most sensitive 
way to study the 2-D distribution of temperature and metallic- 
ity in such a disturbed system as HCG 62. Hi gher -resolution 
spectral maps are presented and discussed in ^5.31 

5. DISCUSSION 
5.1. Interaction of the Radio Plasma with the Hot Gas 
5.1.1. Energetics 

We wish to compare the cavity power - which is a measure 
of the energy injected into the hot gas by the AGN outburst - 
with the gas luminosity inside the cooling radius - which rep- 
resents the energy lost by the hot gas due to X-ray radiation. 
The total power of the two cavities can be estimated directly 
from the X-ray measurement of the pV work done by the jets 
in inflating the cavities, and from measurement of the cavity 
age (e.g., McNamara & Nulsen 2007). By measuring the cav- 
ity enthalpy as H = 4pV (for a relativistic plasma) and the cav- 
ity age as the sound crossing time ts = R/cs = R/ "/kT / ^mn, 
where R is the projected distance from the center of the cavity 
to the group center, 7 = 5/3, and /i = 0.62, we estimate that 
the total cavity power is fcav = 3.8 x 10"^^ ergs"', which is in 
agreement with the results of Birzan et al. (2004) and Rafferty 
et al. (2006). In particular, in deriving the sound speed c,, 
we adopt the temperature values estimated i n the first NE and 
SW sectoral annuli, respectively (see Sect. 14.31 ). The results 
are reported in Table [3] Given the uncertainties in estimating 
the cavity age, for comparison we also adopt the methods ex- 
plained in Birzan et al. (2004) which consider the buoyancy- 
time fbuoy ^ R/ jSC, where g = GM<r/R^ is the grav- 
itational acceleration, V is the volume of the cavity, S is the 
cross-section of the cavity, and C = 0.75 is the drag coeffi- 
cient (Churazov et al. 2001), and the refill-time f,ef ^ 2yjr/g, 
where r is the radius of the cavity. We calculate the gravi- 
tational acceleration at the distance of the cavities from the 
total mass profile derived by Morita et al. (2006), who es- 
timate M<K^8.gkpc ~ 2 X 1O"M0. We find that the age esti- 
mates agree to within a factor of 2, with the buoyancy times 
of ^ 2.0 X 10^ yr lying in between the sound crossing time of 
^ 1.8 X 10^ yr and the refill times of ~ 3.7 x 10' yr. 

We estimate the gas luminosity inside the cooling radius 
('"cool = 33 kpc, Rafferty et al. 2006) as the bolometric X- 
ray luminosity, Lx, derived from a deprojection analysis of 
XMM-Newton spectra similar to that described in ^4.21 We 
measure Lx = 1.5 y. lO'^^ergs"', which is consistent with the 
value of Rafferty et al. (2006). This indicates that the me- 
chanical luminosity of the AGN outburst is large enough to 
balance the radiative losses. In fact, since the 235 MHz radio 
emission extends beyond the cavities, we may argue that there 
is more power in the jet and lobes than one can infer from the 
X-rays. The cavity volumes may be larger as the radio emis- 
sion is less sensitive to projection effects than depressions in 
the X-ray image (as also pointed out by Birzan et al. 2008). 
By assuming that the cavities might be as extended as the ra- 
dio lobes observed at 235 MHz, we estimate that cavity power 
could increase by a factor ~ 2 (see TableO. 

Recent studies have shown that the mechanical cavity 
power usually exceeds the radio luminosity of the bubbles in 
galaxy groups and clusters, with a large scatter from a few 
to a few thousands that seems to show a trend with the radio 
luminosity itself, i.e. the ratio of cavity power to radio lumi- 
nosity increases with decreasing radio luminosities (Birzan et 



TABLE 3 
Cavity properties 



Cavity N 



Cavity S 





(using Vx) 


(using Kadio) 


(using Vx) 


(using Vradio) 


a (kpc) 


5 


10 


4 


10 


h (kpc) 


4.3 


6 


4 


6 


R (kpc) 


8.4 


15 


8.6 


15 


V (cm') 


1.14 X 10*^ 


4.43 X 10" 


7.88 X 10'^'' 


4.43 X 10*' 


p (erg cm"') 


2.5 X 10"" 


1.8 X 10"" 


3.2 X 10"" 


2.2 X 10"" 


pV (erg) 


2.9 X lO^*- 


8.0 X lO'*" 


2.5 X lO^'' 


9.8 X lO"^ 


kT (keV) 


0.75 


0.84 


0.83 


0.95 


fs (kms"') 


440 


465 


463 


495 


h (yr) 


1.9 X 10' 


3.2 X 10' 


1.8 X 10' 


3.0 X 10' 


Pcav(ergs"') 


2.0 X 10*2 


3.2 X lO'*^ 


1.8 X 10^*2 


4.2 X lO'*^ 



Note. — Cavity volumes are calculated assuming spherical or prolate 
ellipsoidal shapes with semimajor axis a and semiminor axis h. In particular. 
X-ray volumes Vx are estimated from Chandra data and are consistent with 
the values of Birzan et al. (2004) and Rafferty et al. (2006), whereas ra- 
dio v olumes V,aciio are estimated from GMRT observations at 235 MHz (see 
j|3.2t . Cavity powers Pcav are calculated assuming ApV of energy per cavity 
and the sound crossing timescale ts- We adopt the temperature and density 
values estimated in the NE and SW sectors from XMM-Newton spectra (see 
top left of Fig. H). 

al. 2004, 2008). Thanks to the new GMRT radio data, we 
can compare the AGN mechanical power with the radio lumi- 
nosity of the source in order to directly estimate its radiative 
efficiency. The monochromatic 235 MHz power of the to- 
tal radio source is P^a = 1.8 x lO^^WHz"'. By assuming a 
spectral index a =1.2 (see ^3.2| i, we measure the total radio 
luminosity over the frequency range of 10 MHz - 10 GHz to 
be 3.0 x 10^^ erg s~' , which is about four orders of magni- 
tude less than the total power of the cavities. Therefore the 
radio source in HCG 62 has a synchrotron radiative efficiency 
as low as ^ 10"'*. This makes HCG 62 the least-efficient sys- 
tem in the Birzan et al. (2008) sample, and is a confirmation 
that the synchrotron radio luminosity is not a reliable gauge of 
the total mechanical power of the AGN outburst, as the radio 
sources can be very poor or time-variable radiators. 

5.1.2. Pressure Balance 

Since the radio source is filling the cavities, we can compare 
directly the radio pressure of the relativistic plasma internal to 
the lobes with the X-ray pressure of the surrounding thermal 
gas. The pressure of the hot gas is measured from the density 
and temperature derived from the X-ray data as = In^kT 
04. 2b . In particular, for the X-ray pressure in the northern 
and southern cavities we assume the values estimated in the 
first bin by performing a deprojection analysis along the NE 
and SW sectors, respectively (see top left of Fig. |5]l. We find: 
/?x,N = 2.5 X 10""ergcm"^, /?x.s = 3.2 X 10""ergcm"-'. These 
values are consistent with the gas pressure at the radius of the 
cavities estimated by Morita et al. (2006). 

The total pressure in a radio lobe is the sum of the magnetic 
pressure, pB, and the total particle pressure, ppart, and can be 
written as 



Pradio - PB+ Ppait - 



part 



B\nil]^ (2) 



(fiV 3 (fiV 87r 3 (fiV 



where k is the ratio of the energy in protons to that in electrons 
{E^), V is the volume of the radio lobe and is the volume 
filling factor Using the expression for Zig given in Pachol- 
czyz (1970), Eq. |2] determines the lobe pressure in terms of 
the magnetic field strength and the factor k/ip, once the vol- 
ume V of the radio lobe is known. This calculation is usually 
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TABLE 4 
Equipartition calculations 



Cavity N Cavity S 

(revised) (standai'd) (revised) (standard) 



5235 (mJy) 


6.0 


6.0 


5.4 


5.4 


5610 (mJy) 


1.8 


1.8 


1.7 


1.7 


610 
"235 


1.3 


1.3 


1.2 


1.2 


P235 (WHz-') 


2.5 X lO^i 


2.5 X 10^' 


2.3 X lO^i 


2.3 X 10^ 


O'miii 


100 


900 


100 


850 




6.7 


3.1 


6.6 


3.3 


Pradio (ergcm"^) 


2.3 X 10-'2 


4.8 X 10-" 


2.3 X 10-'2 


5.5 X 10-' 


Px/Pradio 


10.9 


52.1 


13.9 


58.2 


^bal 


318 


1830 


500 


2450 



nine FR-I radio galaxies in a sample of groups of galaxies, 
Croston et al. (2008) found that the radio lobes are under- 
pressured by a factor up to 70 (with revised equipartition, as- 
suming 7niin = 10 and a = 0.5), and that the pressure imbal- 
ance appears to be linked to the radio-source morphology, i.e. 
'plumed' sources typically have larger pressure deficits than 
sources where the jets are embedded in the lobes ('bridged' 
sources). The authors interpret this result as evidence that 
plumed sources have a higher entrainment rate due to the 
larger fraction of the jet surface which is in direct contact with 
the external medium, leading to an increase in k/cf). Although 
the classification into bridged and plumed morphologies may 
not directly apply to radio sources at the center of cooling 
core systems, typically having amorphous structures, this pic- 
ture is consistent with the results of Dunn et al. (2006) who 
argue that the large pressure imbalance observed in radio bub- 
bles as those of the Perseus cluster is more likely to be due to 
entrainment rather than a relativistic proton population. Un- 
fortunately we cannot derive good constraints on the energy 
in relativistic particles without knowing the break frequency 
in the radio spectrum. However, we note that HCG 62 shows 
a disturbed radio morphology with inner lobes clearly filling 
the well defined X-ray cavities, b ut wi th outer lobes having 
no associated X-ray cavities (see ^3.2l i. Assuming their de- 
tection is not limited by the sensitivity of the current Chandra 
image, this suggests the possibility of mixing between ambi- 
ent gas and radio plasma in the lobes. Therefore the A:bai > 
values that we measure in the lobes might be the results of 
entrainment of thermal gas through the group atmosphere. 

5.2. Shock Fronts 

We wish to investigate in more detail the morphological fea- 
tures presented in 33.1l bv measuring the surface brightness 
and temperature profiles. 

The surface brightness profile extracted along the the SW 
sector within an aperture of 260° - 330° (counting counter- 
clockwise from P.A. = 0° toward the west) is shown in the 
top-left of Fig. |6] By fitting a broken power-law density 
model, at 36 kpc from the center we identify a clear break 
in surface brightess which corresponds to a jump in density of 
1 .65. To unambiguously determine the nature of this front, we 
extracted the temperature profile along the same sector (Fig. 
|6] bottom left) and find that the region immediately interior 
to the front is significantly hotter than the undisturbed region 
just outside of it, with a temperature jump across the front of 
^14%. As discussed below, this is consistent with the inter- 
pretation of the front as a shock. 

We also identify a front to the NE at 20 kpc from the cen- 
ter (see Fig. |6] top right), that appears to have a symmet- 
ric feature to the SW (the surface brightness profile shows an 
edge at 16 kpc from the center, see inner box in Fig. |6] top 
left). We note that the detection of a temperature rise in the 
regions immediately inside of these fronts, which would be 
expected if they are shocks, is complicated by the underly- 
ing rising temperature profile of the global group atmosphere 
(see Fig. |6] bottom). For this reason, despite the lack of evi- 
dence for a temperature jump across the inner fronts, we can- 
not rule out the presence of shocks. In fact, the symmetric 
placement of these fronts and their proximity to the cavities 
(3-5 kpc outside the cavity edge) support their interpretation 
as weak shocks. Alternatively, they could be interpreted as 
cold fronts, caused by the subsonic expansion of the cavities 
displacing cool, low entropy gas outward. However, since we 
cannot place conclusive constraints on these features with the 



Note. — Results of e quipart ition calculations. We assume (p = I, k = I 
and V = Vx estimated in ^5.1.11 A low energy cut-off 7,^1,, = 100 of the 
electron spectrum is assumed to calculate the "revised" equipartition field, 
which implies a spectral frequency range of ~ 280 kHz - 100 GHz. We also 
show for comparison the results of the standard equipartition calculations, 
which adopt instead a fixed emitted frequency band 10 MHz - 100 GHz. 

performed under the widely adopted minimum energy condi- 
tions, in which the relativistic plasma is in equipartition with 
the magnetic field (Beq)- For historical reasons the frequency 
band adopted to calculate the standard equipartition field is 
I'l = 10 MHz -1^2 = 100 GHz, i.e. roughly the frequency range 
observable with radio telescopes. From a physical point of 
view, the adoption of this frequency band in the calculation 
of the minimum energy is equivalent to the assumption that 
only electrons emitting between 10 MHz - 100 GHz, i.e. with 
energy between 7inin oc (yi /Beq)'''^ and 7n,ax oc (yi/Beq)^^^ are 
present in the radio source. This approach neglects the con- 
tribution of the electrons emitting below 10 MHz and, as a 
more serious bias, in radio sources with different Beq selects 
different energy bands of the electron population because the 
energy of the electrons which emit synchrotron radiation at 
a given frequency depends on the magnetic field intensity 
(Brunetti 2002). We therefore adopt a different approach to 
calculate the minimum energy conditions, in which Bgq does 
not depend on the emitted frequency band but directly on the 
low energy cut-off of the electron spectrum. These so-called 
"revised" equipartition conditions select also the contribution 
to the energetics due to the low-energy electrons (Brunetti et 
al. 1997). In our calculations we assume j„^i„ = 100 and the 
observed spectral index (see Tab. |4|i. 

The results of the equipartition calculations are reported in 
Table |4] where for comparison we also show the standard 
equipartition values. We estimate that the X-ray pressure is 
about one order of magnitude higher than the radio pressure, 
as typically found in cavity systems (e.g., Blanton et al. 2001, 
De Young 2006, Croston et al. 2008). We also find that with 
revised equipartition the cavities are closer to pressure bal- 
ance than they are with standard equipartition. On the other 
hand, by assuming the lobes are in pressure equilibrium with 
the ambient gas we can obtain constraints on the particle con- 
tent within the radio lobes (Dunn & Fabian 2004, Birzan et al. 
2008). In particular, we determine the ratio ^bai of the energy 
in protons to that in electrons that is required to achieve pres- 
sure balance under revised equipartition conditions. We find 
^bai ~ 300-500, whereas higher values (^bai ~ few thousands) 
are expected with standard equipartition. These values are in 
the range found by Birzan et al. (2008), who studied the en- 
ergetics and particle content of the lobes of 24 radio galaxies 
in cooling cores, obtaining values of A;bai up to approximately 
4000 (with standard equipartition). Based on an analysis of 



10 




Radius (kpc) radius (kpc) 




100 



radius (kpc) 



radius (kpc) 



Fig. 6. — Top, left: Background-subtracted, exposure-corrected Chandra surface brightness profile extracted along the SW sector from P.A. 260° to 330° in 
the energy range 0.5-2.0 keV. The surface brightness is in counts cm"^ s"', with errors at Icr. Radial error bars show the limits of the bins. The smooth curves 
show fits of the broken power-law density model which give density jumps of 1.47 (dotted), 1.65 (solid), and 1.86 (dashed), corresponding to Mach numbers of 
1.32 (dotted), 1.45 (solid), and 1.62 (dashed). The radius of the shock is 36 kpc. The inner box shows the surface brightness profile of the inner front, which is at 
a radius of 16 kpc. The solid curve shows the fit of the broken power-law density model which gives a density jump of 1 .77. Top, right: Similar to top left panel, 
but for NE sector from P.A. 95° to 160° . The solid curve shows the fit of the broken power-law density model which gives a density jumps of 2.27. The radius of 
the front is 20 kpc. Bottom, left: Chandra (black circles) and XMM-Newton (green triangles) temperature profiles measured along the SW sector (same sector as 
in top left panel). The dashed line indicates the position of the shock (rshock = 36 kpc), while the dotted line indicate the position of the inner front (rfront.sw = 16 
kpc). Bottom, right: Chandra temperature profile measured along the the NE sector (same sector as in top right panel). The dotted line indicates the position of 
the front (ffront.NE = 20 kpc). Error bars on temperature values are at la. 



present data, we will not discuss them further in the paper 
and in the following we will focus only on the outer shock 
detected to the SW (simply referred as to 'the shock'). 

A hydrodynamic model for the shock was made by initi- 
ating an explosion at the center of a hydrostatic, isothermal 
atmosphere with a power-law density profile. The power-law 
index for the density for the unshocked gas was determined 
from the broken power-law fit to the surface brightness pro- 
file. Our data are consistent with the presence of a shock 
having a Mach number = 1.45 and an age >2.7 x 10^ yr 
Compared to a model with continuous energy injection, the 
point explosion produces a stronger, i.e. faster, shock at early 
times so its age estimate is a lower limit. For such a Mach 



number our model predicts the emission measure weighted 
temperature to rise by ^\5% across the front, which is in 
agreement with our measurements. As it appears evident 
from the comparison with the temperature profile obtained 
with XMM-Newton along the same sector (overlaid in green 
in Fig. |6] bottom left), the detection of the temperature jump 
is possible only with the superb spatial resolution of Chan- 
dra. We estimate an energy ^ 3.4 x 10^** erg, and a power 
^ 4.0 X 10"^^ ergs"' which is about one order of magnitude 
higher than the cavity power The unusually large difference 
between the shock power and the cavity power is essentially 
due to the fact that the shock is much further away than the 
cavities from the group center (./^shock = 36 kpc vs. /?cav ^8.5 
kpc). Such a difference may be reduced by a factor ^ 2 if we 
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Fig. 7. — Left: Chandra temperature map with GMRT 235 MHz radio emission overlaid in black contours. Superposed in white are the two elliptical regions 
indicating the cavities (thin), and two arcs indicating the SW shock front (thick) and the NE front (thick dashed). Right: XMM temperature map obtained by 
means of the hardness ratio of difi'erent energy bands. Overlaid radio contours and superposed regions are the same as in left panel. Owing to the methods 
adopted to build such maps, the pixel values are not independent (see text for details). 



adopt the cavity volume estimated from the radio data instead 
of that estimated from X-rays. 

From a physical point of view, the difference in location 
and energetics between the shock and cavity system can be 
explained by two possibilities: 1) the shock and cavities are 
created by the same, violent AGN outburst; 2) the shock and 
cavities originate from different, multiple episodes of gentle 
AGN outbursts. 1) In the first case, the jet is initially super- 
sonic and inflates the cavity violently, with much of the energy 
driving a shock. Once the jet slows the cavity rises buoyantly 
and the shock detaches from its tip and continues to propa- 
gate in the group atmosphere at a velocity A^c, . By assuming 
that the shock was never any weaker than it is now and that 
the group atmosphere has an average temperature of 1 ke V, 
we estimate the time for the shock to travel from the AGN to 
its current location to be <4.8 x 10' yr, which is consistent 
with the shock age estimated above. The shock age is also 
consistent with the refill timescale of the cavities, but slightly 
longer than the buoyancy timescale (see ^5. 1.1b . The buoy- 
ancy timescale can be considered as an upper limit to the age 
of the outburst that produced the cavities, since it assumes 
the cavities were inflated at the group core and have risen 
slowly. If they formed close to their current location their 
age would be even shorter, and the conflict with the shock 
age more serious. 2) In the second case, it is possible that 
the shock originated during a previous AGN outburst which 
produced the outer radio lobes and perhaps an associated pair 
of cavities, though these are not detected in the existing data. 
In this scenario, the inner cavities are thought to be evacuated 
by the inner radio lobes during a more recent outburst, thus 
providing a natural explanation for the marginal difference 
between the buoyancy timescale of the cavities and the age 
of the shock. The expected older outer cavities might be de- 
tected with deeper exposures in a region close to the observed 
shock front. 

However, with the current radio data we are not able to 



determine the relationship of the outer radio lobes to the in- 
ner ones. Therefore, considering also the uncertainties in the 
timescales estimated above, we cannot determine whether the 
observed radio/X-ray morphology is the result of one single 
AGN outburst, or multiple episodes. In either case, the posi- 
tion of the shock outside the S radio lobe makes plausible the 
interpretation of the shock as being directly driven by the lobe 
expansion triggered by an AGN outburst. Averaging over the 
past few 10' yr, the cavity power alone thus provides a lower 
limit to the true total mechanical power of the AGN. Inclusion 
of the energy in the shock provides a more complete estimate. 
The shock plus cavity power of ~ 4.4 x 10'*"' ergs"' is much 
higher than the co oling luminosity {Lx = 1.5 x 10'*^ ergs"', es- 
timated in ^5. 1.1b . If the power level of 4.4 x 10"*^ ergs"' is 
sustained, since it exceeds the power radiated by cooling core, 
either the core is currently being heated, or most of the power 
must be deposited in gas beyond the cooling core. The fact 
that the shock front is located at the outer edge of the cooling 
region is consistent with the latter 

On the other hand, the apparent absence of an outer shock 
front to the NE direction, which would be expected if the 
shock was created by the two outer radio lobes, requires some 
asymmetry. It is also possible that the shock strength varies 
with the position angle. In this case the shock front may be in 
fact intrinsically symmetric and surrounding the whole central 
region of the group, but is detected only to the SW because of 
some difference in gas properties (e.g., higher density). We 
investigate further these possibilities also by means of a tem- 
perature map, which shows a lack of very cool (<0.9 keV) 
material on the SW side of the core (see ^5.3.11 ) that might be 
explained by the shock. 

5.3. Spectral Maps 

Since no obvious feature associated w ith th e cavities and 
radio lobes appear in the wedge maps 04.31) . we produced 
spectral maps with higher spatial resolution in order to investi- 
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gate the 2-D temperature 05.3.1I ) and iron abundance ( ^5.3.2l i 
distribution independently of the choice of the sectors used to 
extract the spectra. Such spectral maps are suitable for identi- 
fying potentially interesting structures. 

Chandra spectral maps were produced using a method sim- 
ilar to that described in O' Sullivan et al. (2005). The ACIS- 
S3 field-of-view was divided into a grid of square map pixels 
with side length 10 physical pixels ('^5"). For each of these 
map pixels a spectrum was extracted from a circular region 
centered on the pixel, with a radius chosen so as to obtain 
800 net counts in the 0.3-5.0 keV band. Radii were allowed 
to vary between 10 and 100 physical pixels. Spectra and re- 
sponses were created in an identical fashion to that used in 
our main spectral analysis and fitted with an absorbed apec 
model. Point sources were excluded from both spectra and 
the calculation of the spectral extraction region size and en- 
ergies outside the 0.3-5.0 keV band were excluded. 90% un- 
certainties on parameters were estimated and any pixel with a 
temperature uncertainty greater than 15% was excluded. The 
colors of the remaining map pixels indicate the best fitting 
temperature, or abundance. 

5.3.1. Temperature 

In Fig. |7](left) we show the temperature map obtained from 
Chandra spectra. We also show in Fig. |7](right) a temperature 
image of the group central region built from XMM-Newton 
colors. Specifically, we extracted mosaiced MOS images in 
four different energy bands (0.5-1.0 keV, 1.0-2.0 keV, 2.0-4.5 
keV and 4.5-8.0 keV), subtract the background and divide the 
resulting images by the exposure maps. A temperature in each 
square map pixel, chosen to have side length of 5 physical 
pixels (= 5.5"), is then obtained by fitting a thermal plasma 
model to the count rate estimated in the above energy bands 
in a surrounding area with a fixed radius of 15 physical pixels 
(= 16.5"). 

The two temperature maps are consistent with each other 
and with previous results (Morita et al. 2006; Sanders et al. 
2009), showing the presence of the cool core, but also the 
presence of cool gas in the region of the northern cavity. Ow- 
ing to emission weighting in the large map extraction regions, 
this cold feature could be related to the presence of the cavity 
rims, which are indeed typically observed to be cold (McNa- 
mara & Nulsen 2007, and references therein). On the other 
hand, the lack of a comparable cold feature in the region of 
the southern cavity is notable. One possible explanation is that 
there was initially the presence of cold gas also to the south, 
but that this has then been heated by the passage of the shock. 
This argument could be taken as an indication of an asym- 
metric shock, or at least of an asymmetric shock strength, as 
discussed in ^5.21 

5.3.2. Iron abundance 

In Fig. |8]we show the iron abundance map obtained from 
Chandra spectra. The iron abundance distribution appears to 
be inhomogeneous, asymmetric, and lacking the usual cen- 
tral peak. We have reason to believe that the structures seen 
in the abundance map are real, as they are also observed in 
the maps of 90% upper and lower bound on abundance (not 
shown here). However, the uncertainties on the value of in- 
dividual map pixels is large (30-50%) and the variable sizing 
of spectral extraction regions produces an effect analogous to 
adaptive smoothing. It is therefore necessary to extract sep- 
arate spectra to determine the true abundance of individual 
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Fig. 8. — Chandra iron abundance map with GMRT 235 MHz radio emis- 
sion overlaid in black contours. Superposed in white are the two elliptical re- 
gions indicating the cavities (thin), and the two arcs indicating the SW shock 
front (thick) and the NE front (thick dashed). 

regions and the significance of any abundance differences. In 
particular, to confirm the presence of the arc-like region of 
enriched material visible at about 2' to the SW of the group 
center, first discovered by Gu et al. (2007), we extract spec- 
tra from two elliptical regions centered at 2' to the NE and 
SW of the group center along the direction of the inner radio 
lobes. The semimajor axis of both ellipses is orthogonal to 
the direction of the inner radio lobes and 1' in extent, with the 
semiminor axis being half of it. We measure abundances of 
0.12!{!:J]2 solar and 0.95![;:?^ solar (la errors) in the NE and 
SW regions, respectively. This indicates that the SW arc is 
more abundant than the NE dip at Sa significance. 

Gu et al. (2007) discuss several possibilities about the ori- 
gin of the high-abundance gas, including entrainment during 
the cavity expansion triggered by the AGN activity, and ram 
pressure stripping due to tidal interaction of two central galax- 
ies which have experienced a recent (^ 10** yr) minor merger, 
as indicated by optical studies (Spavone et al. 2006). Unfor- 
tunately, we do not have enough information in the current 
data to disentangle this problem and the origin of the high- 
abundance arc-like region is still an open question. How- 
ever, we stress here that the enriched material has a loca- 
tion and shape consistent with that of the newly discovered 
shock. Recent papers by Kaastra et al. (2009) and Prokhorov 
(2009) claim that in some cases, for instance near interfaces 
of hot and cold gas and near shocks, the usual approxima- 
tion of Maxwellian electron distributions adopted for calcu- 
lating thermal X-ray spectra is no longer valid. In such sit- 
uations, the presence of non-thermal electrons modifies the 
line emissivities hence the X-ray spectrum, thus affecting the 
measurements of metal abundances. In particular, the best-fit 
iron abundance for the isothermal model is about 30% higher 
than the actual abundance (Kaastra et al. 2009). This could 
explain the high abundance arc that we measure. A similar 
effect of apparent iron overabundance may also be produced 
by the presence of non-equilibrium ionization states below the 
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shock (Prokhorov 2010). These interpretations of the arc-like 
feature observed in the abundance map of HCG 62 are con- 
sistent with the presence of the detected shock front. 

6. SUMMARY 

We have analyzed new low-frequency GMRT radio data and 
existing Chandra and XMM-Newton X-ray data of the com- 
pact group HCG 62 in order to study the properties of the 
cavity system and the interplay with the central radio source. 
Our investigation demonstrates the power of a combined X- 
ray/radio approach to the feedback problem, and particularly 
the benefits of extending radio studies of AGNs to low fre- 
quencies where less energetic, older electron populations are 
visible. We summarize our main results below. 



ergy in protons to that in electrons is ^-^300-500, which 
may originate from thermal gas entrainment in the lobes 
through the group atmosphere. 

4. We discover a Mach ^1.5 shock front at 36 kpc from 
the center toward the SW, which has not been reported 
in literature. The total energy in cavities and shock is 
-3.6 X 10^'^ erg. 

5. We confirm the presence of an arc-like region of en- 
riched gas to the SW whose apparent high metallicity 
could be due to a non-Maxwellian electron distribu- 
tion or non-equilibrium ionization states near the shock 
front. 



1 . In contrast with the early classification of the HCG 62 
cavity system as "radio ghost", the new GMRT obser- 
vations detect low-frequency radio emission in the cav- 
ities. 

2. The radio source has a total radio luminosity of 3.0 x 
10"^^ erg s~' and a very low radiative synchrotron effi- 
ciency of — 10""^. 

3. We estimate that the X-ray pressure of the cavities is 
about one order of magnitude higher than the "revised" 
equipartition pressure of the radio lobes. By assum- 
ing pressure balance, we find that the ratio of the en- 
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